Anharmonic magnetic deformation of self-assembled molecular nanocapsules 
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High magnetic fields were used to deform spherical nanocapsules, self-assembled from bola- 
amphiphilic sexithiophene molecules. At low fields the deformation - measured through linear 
birefringence - scales quadratically with the capsule radius and with the magnetic field strength. 
These data confirm a long standing theoretical prediction (W. Helfrich, Phys. Lett. 43A, 409 
(1973)), and permits the determination of the bending rigidity of the capsules as (2.6±0.8) x 10~ 21 
J. At high fields, an enhanced rigidity is found which cannot be explained within the Helfrich model. 
We propose a complete form of the free energy functional that accounts for this behaviour, and allows 
discussion of the formation and stability of nanocapsules in solution. 

PACS numbers: 62.25.+g, 81.16.Fg, 82.70.Uv, 83.60.Np 



The spontaneous formation of hollow molecular 
spheres (vesicles) is a prototype of self-assembly of or- 
ganic molecules into well-defined architectures. Vesicu- 
lar self-assembly is a quite general phenomenon [l| ob- 
served in lipid molecules (the building blocks of cell 
membranes 2|), polymers, dendrimers, and 7r-conjugated 
oligomers [||. This large variety makes vesicles inter- 
esting for applications like tailor-made nanocapsules for 
drug delivery [l|, opto-electronics [J], and containers for 
chemical reactions [5j . They also provide ideal systems to 
develop physical models for self assembly. Such a model 
requires a quantitative understanding of the weak non- 
covalent interactions involved in the self-assembly process 
on the basis of which vesicle formation and stability can 
be explained. To this end, we have measured magnetic 
deformation of novel bola-amphiphilic @ sexithiophene- 
based capsules 0, Q and we have developed a theory 
to describe this effect. Comparison between theory and 
experiments allows the determination of important quan- 
tities like the bending rigidity of the capsules and con- 
tribute to an understanding of their stability. 

The determination of the elasticity of vesicles through 
magnetic deformation was proposed three decades ago 
by Helfrich 0,[T3, but it was never demonstrated exper- 
imentally. A magnetic field is a thermodynamic param- 
eter that exerts well-defined magnetic forces to deform 
all capsules in solution, in a controlled and contact-free 
way. Furthermore, the field-induced deformation can be 
quantified via a linear birefringence measurement. The 
results are, therefore, easier to compare with theory than 
those of other deformation methods, such as using cap- 
illary forces [TTI . [l2l [p3 | , pressure 0, HU , electric fields 



17 1 , mechanical forces by an atomic force microscope 
or optical tweezers [19(. 
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Theoretically we use a fourth order Landau-Ginzburg 
expression of the capsule free energy that describes quan- 
titatively the measured elastic deformation up to high 



magnetic fields, far beyond the quadratic dependence 
predicted by Helfrich at low fields [10]. From the data we 
obtain the free energy as a function of temperature that 
explains the overall stability of nanocapsules and predicts 
their dissolution at high temperature. The quantitative 
determination of the free energy can be used to develop 
and validate models for the non-covalent interactions re- 
sponsible for the self-assembly process. 

Experimentally we investigate sexithiophene (6T) 
(2,5"" '-(R-2-methyl-3,6,9,12,15-pentaoxahexadecyl es- 
ter), Fig. la) which is a 7r-conjugated oligomer with 
a rigid apolar sexithiophene block substituted at both 
ends by polar ethylene oxide chains. 6T forms hollow 
spherical self-assemblies in a 2-propanol solution (Fig. 
lb, details can be found in Q). Dynamic light scattering 
(DLS) reveals a narrow distribution of capsule sizes, 
with an average radius R(T) (±10%) that increases with 
temperature T from 56 nm at 20 °C to 126 nm at 60 °C 
(Table UJ). At 70 °C the spherical shape is lost, while 
at 80° C the scattering intensity is zero because the 6T 
molecules are molecularly dissolved j§[ . 

6T molecules have a strongly anisotropic diamagnetic 
susceptibility x-, which means that the sexithiophene core 
tends to orient along the field [2(j. Consequently, cap- 
sules are deformed from a sphere to an oblate ellipsoid, 
where more molecules are parallel (top and bottom of 
the capsule) than perpendicular to the field (around the 
equator, Fig. lc), reducing the magnetic energy at the 
expense of elastic energy. This magnetic deformation is 
probed by measuring the field induced optical birefrin- 
gence An in a resistive 20 T magnet using a laser (HcNc, 
632.85 nm) modulation technique. 

Fig. Id shows An of a 1 g/1 6T/2-propanol solution 
at different T. These curves were recorded during slow 
sweeps of the magnetic field B at constant T (±0.1 °C) 
and are fully reproducible for up-, down-, and consecu- 
tive sweeps. In our set-up a negative An corresponds to 
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FIG. 1: (color online) a) Chemical structure of a sexithio- 
phene molecule (6T). Schematic representation of a spherical 
capsule with radius R (b) and magnetically deformed capsule 
with semi-axes a and c (c) d) Temperature dependence of the 
magnetic field (B) induced birefringence of a 1 g/1 6T solution 
in 2-propanol. Inset: normalized birefringence Eq.(|3]) divided 
by R 2 and plotted versus B 2 . The dashed line corresponds to 
a fit with constant bending rigidity k. 



6T molecules aligning with their long axis parallel to the 
magnetic field direction, as expected for an oblate ellip- 
soid (Fig. lc). The absolute value of An increases with 
T up to 60 °C. At 70 °C it is considerably reduced with a 
different field dependence, and An is zero at 80 °C, when 
the capsules are dissolved. At low fields An scales with 
the square of the capsule radius and increases as B 2 with 
approximately the same slope at all temperatures (in- 
set of Fig. Id). These experimental facts prove that An 
originates from a reversible effect on stable nanocapsules 
below their melting temperature (20 °C-^60 °C). 

For small deformations (low fields) the normalized bire- 
fringence An/An max (An max = —6.5 • 10~ 5 [2lj]) varies 
linearly with deformation 0, [l(| , but a more accurate 
relation is needed for high deformations (high fields) . For 
small capsule concentrations p c the refractive index n is 
related to the polarizability of individual capsules a by 
n 2 = n,Q + 4irp c a, where no is the refractive index of the 
solvent. The dipole moment of the 6T molecules p is con- 
nected with the effective electric field E' in the molecular 
layer by p = a± (E'— e (eE'))+a||e (eE'), where a± and 
«H are the polarizabilities of 6T molecules with respect to 
e, the unit vector normal to the capsule at a given point 
of its surface. Applying standard boundary conditions 
on the Maxwell equations we find the birefringence of a 



uniaxially deformed capsule: 

An = np c p m ( ay — - a 



dS(l-3e 2 ) 



(1) 



where e and e n are the dielectric constants of the solvent 
and of the 6T in the direction along e, respectively, p m 
is the surface density of molecules in the layer, the inte- 
gration is over the surface S of the capsule and z labels 
the direction of magnetic field B. 

To calculate An we assume that the shape of the de- 
formed capsules can be described by superspheroids, pa- 
rameterized by : 



x(u,v) = a ■ sin" v cos it, 
y(u,v) — a ■ sin' v sin u, 
z{u, v) — c ■ cos 1 v, (a > c), 



(2) 



with a is the radius in the circular xy plane and c the 
extension in the z direction (Fig. [Tk). Superspheroids 
are more general than the spheroids (I = 1) that are 
commonly used for small deformations [l(| , and allow 
the capsules to flatten (I < 1). The birefringence now 
reads: 
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(3) 



The inset of Fig. [2] shows An/An max as function of de- 
formation (a — c)/R for (supcr)spheroids. We keep the 
surface area of the capsule constant, assuming a constant 
number of 6T molecules that form a permeable mem- 
brane. Remarkably, both curves are rather similar, in 
particular in the low deformation limit, where we recover 
the linear relationship between An and deformation by 
expanding Eq. © as HdSl: 



A? 



An n 



R ' 



(4) 



which holds as long as (a — c)/R < 0.3 (inset Fig. [2]). 

We determine the equilibrium values of a, c and I at 
each B by minimizing the total free energy Ttot , which is 
the sum of the elastic J- e \ and magnetic ^"mag free energy. 
Simultaneously, the parameters in J- e \ are determined by 
a fit to the experimental An using relation ([3]). We first 
consider the well known form of the elastic term proposed 
by Helfrich [l(|, leading to: 



•^tot — Pel + J mag + ^surf = — I dS H 2 — 



7 dS-AwR' 



(5) 



where H = (n\ + K2)/2 is the mean curvature, K\, k,2 are 
the principal curvatures, D — 6.4 nm (length of long axis 
of 6T) is the membrane thickness, and Ax = X-l ~ X||> 
the difference in magnetic susce ptib ility along the short 
and long axes of 6T respectively [231 ] . The term ^surf im- 
poses, for large 7 (24|, the constraint of constant capsule 
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FIG. 2: (color online) Normalized birefringence of magneti- 
cally deformed 6T capsules. Solid lines: experimental data 
for different temperatures. Dashed lines: calculated birefrin- 
gence in the Helfrich model, fitting the low-field behaviour us- 
ing a constant k for all T. Inset: calculated relation between 
birefringence and deformation for superspheroids (solid) and 
spheroids (dashed) compared to An/An max ~ (a — c)/R 
(dash-dotted line). 



TABLE I: Structural properties of 6T nanocapsules in 2- 
propanol at different temperatures. Capsule radius R, as de- 
termined by dynamic light scattering (exp.) and obtained by 
fitting the magnetic birefringence (fit). Equilibrium parame- 
ters a, c and I of the deformed capsules at 20 T. 



Temperature 


Radius R 


(nm) 


Superspheroid(20 T) 


T (°C) 


Exp. 


Fit 


a/R 


c/R I 


20 


56 


58 


1.19 


0.59 0.90 


30 


68 


62 


1.20 


0.56 0.89 


40 


75 


75 


1.23 


0.50 0.88 


50 


88 


99 


1.26 


0.42 0.85 


60 


126 


110 


1.28 


0.39 0.84 



surface. The bending rigidity k is the only fit parameter 
to obtain the best description of the experimental bire- 
fringence (Fig. [2]). This approach can be used because 
the capsule size is large compared to the intermolecular 
distances, which means that the capsule shape is charac- 
terized by its curvatures. 

Considering first the limit of small deformations we can 
expand Eq.© as: 



R 



B 



2 R 2 D^X 
3k/io 



(6) 



Combination of Eqs. (J3J and ([6]) shows that the low 
field birefringence is inversely proportional to the bend- 
ing rigidity k, and scales quadratically with magnetic 
field and capsule radius (the Helfrich result). Our ex- 
perimental results (Fig. [2l inset of FigJTj) indeed exhibit 
this scaling up to B ~ 1 T, allowing to determine k 
as (2.6±0.8) x 10~ 21 J, independent of temperature. Our 
data, therefore, permits the experimental determination 



of the bending rigidity of this type of capsules. 

The Helfrich free energy thus quantitatively describes 
the birefringence data in the low deformation regime us- 
ing a constant bending rigidity. However, this model 
clearly overestimates the experimental birefringence at 
high fields (Fig. [J), which indicates an increased bending 
rigidity at higher deformations. Therefore we extend the 
expression of the free energy by including all symmetry 
allowed terms up to fourth order (25|: 

T eX = J dS[A'K + AH 2 + c x H A + c 2 H 2 K + c 3 K 2 ], (7) 

where K = k±K2 is the Gaussian curvature. Note that 
terms in H K and H 3 are not allowed for symmetric bola- 
amphiphiles [25[ . Since at B = the capsules are spher- 
ical, with H 2 — K = 0, it is convenient to rewrite J- c \ as 

jF el = J dS[(A+A r )K+A(H 2 -K) + (c 1 +c 2 +c 3 )H 4 - 

- (c 2 + 2c 3 )H 2 (H 2 -K) + c 3 (H 2 - K) 2 } . (8) 

If we compare the elastic energies Eq.© for a flat layer 
and for a set of spheres with radius R with the same 
total area, and neglecting the edge energy due to the 
interaction of the solvent with the hydrophobic cores, 
we find that A + A' < and c\ + c 2 + c 3 > must 
hold to favour the formation of spheres. In this case, the 
minimum of Eq.© occurs for a sphere of radius: 



R 2 



1 

H 2 



, Ci + c 2 + c 3 
' A + A' ' 



(9) 



Since the term JdSA'K does not depend on deformation, 
the magnetic deformation data cannot provide directly 
information about A'. In fact, according to the Gauss- 
Bonnet theorem, JdSK = Ait for any surface topologi- 
cally equivalent to the sphere (like the superspheroid con- 
sidered here). However, Eq.@ can be used a posteriori 
to determine A! ', once the other parameters are found by 
a best fit to the experiments. This approach is very pow- 
erful, because it experimentally determines the total free 
energy of nanocapsules, which is very difficult to obtain 
otherwise, and which can be used to discuss the overall 
stability of capsules. 

Within this description the low field bending constant 
k can be calculated analytically for spheroids by expand- 
ing Eq.© for small deformations: 



2A 



28ci + 22c 2 + 16c 3 
3^ ' 



(10) 



which connects our approach to that of Helfrich. More- 
over, since k has to be positive to guarantee stability 
against deformations, this equation poses an extra con- 
straint on the choice of parameters. In the spirit of Lan- 
dau theory, we assume that only A depends on T and 
that the coefficients Ci are temperature independent. In 
fact, by using the experimental value of k we eliminate 
the parameter A from the fitting procedure, by which c; 
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FIG. 3: (color online) a) Fit (dashed lines) of the experimen- 
tal birefringence (solid lines) at different temperatures, using 
superspheroidal capsules with deformation energy T B \ given 
by equation (JSj . b) Equilibrium shapes of the deformed cap- 
sules at 20 T. 



and R(T) arc determined via fitting to the experiments 
at the different temperatures (Fig. [S^)- 

For all temperatures the overall behavior of An is 
nicely reproduced, including the weaker field dependence 
at higher fields, indicating anharmonic deformation char- 
acterized by a field enhanced bending rigidity. The fitted 
value of R(T) (Table 0} is well within the error margins 
of the experimental DLS values. Table U also shows the 
parameters a, c and I at 20 T, which correspond to sub- 
stantially deformed capsules, as illustrated in Fig. [5Jd. 
At the highest fields deformations (a — c)/R as large 
as 0.9 are found, which agrees very well with the ob- 
served deformation of 6T capsules that are fixed in a 
compatible organogel and imaged by scanning electron 



microscopy |8j. Not surprisingly, we find that with in- 
creasing deformation the capsule shape becomes flatter 
(decreasing I, see Table HJ, maximizing the number of 
molecules aligned along the field. 

The best fits in Fig. [3^ are obtained using c\ = 
1.148, c 2 = -1.336, c 3 = 0.358 (10~ 35 Jm 2 ), and the In- 
dependent A and A + A' are shown in the inset. In- 
deed, these values make the capsules stable at B = 
(A + A' < 0, ci + c 2 + c 3 > 0) for T < 60°C. How- 
ever, A + A' approaches when T — > 60°C, suggesting 
an instability of the capsules at higher T, which is indeed 
compatible with the experiments that show that the cap- 
sules melt above 60°C. It is this vicinity of the system 
to an instability, which makes fourth order terms com- 
parable in magnitude to second order ones and make the 
Helfrich model inadequate for this situation. 

The parameters c±,C2,Cs, A(T) and A'{T) fully de- 
scribe the free energy of the nanocapsules. Their actual 
values are intricately related to the intermolecular inter- 
actions within the capsule/solvent system, such as ir-ir in- 
teractions between neighboring 6T cores and interactions 
between 6T tails and solvent molecules [Hi, 26] . Determi- 
nation of these parameters for a physical nanostructure 
therefore provides a stringent test case, on the basis of 
which microscopic intermolecular interaction models can 
be developed. Such a quantitative description of non- 
covalent interactions within molecular assemblies is an 
important step towards the detailed understanding of the 
rules of self-assembly. 
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